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1. Introduction 

In MEG and EEG source localization, an ill-posed 
inverse problem has to be solved. This implies that 
in order to obtain a meaningful solution, the data 
alone do not suffice: a source model describing 
properties of the expected sources must be supplied 
as well. 

Two commonly used families of source models are 
the dipole models and the distributed source models. 

1.1 Dipole models 

For dipole models, a relatively low number of 

parameters are adapted to the data. These are the 
dipole locations and moments. The choice of 

parameters is motivated by external knowledge 
about the data to be analyzed, such as the number 
and temporal characteristics of expected sources. 
Due to the low number of parameters (unknowns), 
the solution is often the set of values that explain as 
much of the measured data as possible, although 
regularized solutions are used as well [1]. 

Dipole models have two main disadvantages: 

1. A dipole is an abstract concept that represents 
an extended patch of activated cortex. If the 
patch of activated cortex is curved, the dipole 
may not even be located on that patch but rather 
inside the enclosed volume. 

2. The number of dipoles and their characteristics 
are not always well known a priori. 

1.2 Distributed source models 

Distributed source models (also called current 
density models) assume the simultaneous activity of 
many sources, which are distributed throughout the 
volume-of-interest. Thus, they can theoretically 
resolve both disadvantages of the dipole fit 
methods: 

1. By distributing sources within the cortex, 
cortical activity can be modeled without 
abstraction. 

2. Due to the multitude of sources, the number of 
active configuration is an outcome rather than a 
parameter of the method. 
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Figure 1: a) Cortical source locations; b) dipole 
point sources oriented perpendicular to the cortical 
surface. 


However, other problems arise: 

• The model must describe properties of the 
expected sources that allow to pick a favored 
solution. Solutions of minimum norm or 
maximum smoothness are well-established here. 

• Depth-weighting (lead field normalization) 
techniques have to be used in order to achieve a 
constant sensitivity profile over as much of the 
volume-of-interest as possible. 

• The best-fitting solution will normally over-fit 
the data. This means, that it is normally possible 
to explain more of the data than its signal 
portion only, normally resulting in chaotic 
source distributions. Regularization techniques 
establish equilibrium between closeness to the 
model and goodness-of-fit. 

It can be seen, that depending on the source model, 

there are several trade-offs: 

• While minimum L2 norm methods are fast with 
respect to Lp norms with p < 2, they suffer from 
an artificial smearing of the reconstructed 
source distribution. 

• On the other hand, with p close to 1 and 
vanishing smearing artifacts, the method 
becomes very sensitive to noise in the data [2]. 

• Methods explicitly imposing smoothness 
constraints such as the LORETA method [3] are 
computationally rather slow and perform 
unsatisfying, if smoothness is measured along 
the cortical surface [4]. 

An extended source model might resolve some of 

these trade-offs. 




1.3 Extended sources 

Traditionally, current density models use point 
sources (elementary dipoles) as the basis of their 
computations. In this paper, an extended source 
model is proposed instead. In the extended source 
model, overlapping patches of activity rather than 
point sources enter the - otherwise unchanged - 
current density algorithms. 

2. Methods 

2.1 Point sources 

Standard distributed source models assume the 
simultaneous activity of several thousand point 
sources j. A data and a model term are minimized 
simultaneously. For least squares methods, we get 

j = arg min ||C(m - Lj)|| 2 + X ||Wj|| 2 

= ML r C r (CLML r C r + Al) 1 Cm 

where 

M = W r W 

If sources are known to be extended spatially or 
temporally, this can be modeled by a coupling of 
neighbored activity via a non-diagonal matrix W. 
This approach leads to long computation times due 
to the inversion of W T W. 

We propose a method that models spatial extension 
using extended sources. 
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Figure 2: a) point sources; b) extended sources. 

2.2 Extended Sources 

Extended sources (figure 2) are implemented by 
folding resp. mollification [5]. The mollification 
operator K encodes the relation between extended 
sources and point sources (dipoles) and can thus 
also be seen as a set of point-spread functions. 

A vector u of extended sources is related to j by 
j = Ku. The matrix K encodes the shape of the 
extended sources. For least squares methods, we get 

u = arg min |C(m - LKu)|| 2 + X ||Wu|| 2 

= MK r L r C r (CLKMK r L r C r + X \)'Cm 

Here, W is a diagonal matrix used for depth and 
location weighting and W T W can be inverted on the 
fly. The term LKu ensures that the lead field matrix 
can still be set up for point sources. After u has been 


computed, j is derived and displayed. The matrix K 
is stored and applied in a format that is optimized 
for its sparseness. Possible designs for K include 
gaussian, linear, or boxcar-shaped falloff. An 
extension to the temporal domain is straightforward. 
Then, K can model temporal smoothness or a set of 
activation functions. In all these cases, the minimum 
norm constraint is imposed upon the extended 
sources. 

In this paper, gaussian falloff has been realized by 
successive application of kernels. 

3. Application to epileptic spike data 

For the same patient, different epileptic spikes have 
been recorded using 28-channel EEG and 37- 
channel MEG. Due to the small number of sensors 
and no averaging, these data are a challenge for 
every source reconstruction algorithm. 

3.1 EEG 

Figures 3 and 4 show the sensor setup, the measured 
data, and the volume conductor model used. The 
automatically generated [6] BEM (boundary 
element method) volume conductor model consists 
of three shells with 1275 nodes (skin), 997 nodes 
(skull), and 1326 nodes (brain). 



Figure 3: EEG data (left), contour maps from -25ms 
to 30ms (top right) and volume conductor model 
used (bottom right). 




Figure 4: EEG sensor setup 



















































The segmented and triangulated cortical surface 
consists of 18107 locations (36210 triangles). 

It can also be seen, that due to the small number of 
electrodes, the temporal lobe area is not sampled 
very well. 

For the 5ms time point, the following source 
reconstruction methods have been performed: 

1. Single dipole fit. 

2. Minimum Norm Least Squares (MNLS) 

3. Minimum Lp norm with p= 1.2 

4. Minimum LI norm (linear programming) 

5. LORETA 

6. Min. LI norm with maximum current density. 
Methods 2 to 5 were performed with and without 
using extended sources. Extended sources had a 
diameter of 30mm. For all current density methods, 
the goodness-of-fit was automatically adjusted using 
the % 2 -criterion [5] to be 1/SNR = 7.5%. The 
goodness-of-fit of the single dipole fit was 14.5%. 
Results are shown in figures 5 and 6. 

3.2 MEG 

Figures 7 and 8 show the sensor setup, the measured 
data, and the volume conductor model used. The 
BEM model consists of three shells with 1280 nodes 
(skin), 1010 nodes (skull), and 1330 nodes (brain). 
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Figure 5: EEG reconstructions based upon point 
sources (left column) and extended sources (right 
column). 


Figure 7: MEG data (left), contour maps from 
-52ms to 14ms (top right) and volume conductor 
model used (bottom right). 



Figure 8: MEG sensor setup 



Dipole Fit Minimum LI Norm, max curr. 


Figure 6: EEG dipole fit result (left) and LI norm 
reconstruction with a maximum current density 
constraint (right). 

























































The segmented and triangulated cortical surface 
consists of 18415 locations (36826 triangles). 

The MEG recording was performed for the same 
patient, but a different spike. For the -35ms time 
point of that spike, the same source reconstruction 
methods as in the EEG case have been performed. 
For all current density methods, the goodness-of-fit 
was automatically adjusted using the % 2 -criterion [5] 
to be 1/SNR = 5.5%. The goodness-of-fit of the 
single dipole fit was 11.9%. Results are shown in 
figures 9 and 10. 

4. Results and Discussion 

The extended source model uses source patches 
rather than dipoles as the basis of all computations. 

It reconstructs extended sources in little more than 
the time needed to perform a standard solution. 
Source patches include cortical locations of varying 
orientations, representing a net current density. 
Because these patches overlap, it is still possible to 
localize patches with a resolution in the order of the 
cortical discretization, i.e. a few mm. 

The extended character of the results obtained with 
this method is a prior, not an artifact produced by 
the underlying norm, as with minimum norm least 
squares (MNLS). 

When the method is used with the LI norm, which 
normally produces few point sources, few extended 
sources are computed. The minimum LI norm 
method seems to be an appropriate framework for 
extended sources, as this model does not produce 
smooth source configurations and the method is 
very sensitive. 
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Figure 9: MEG reconstructions based upon point 
sources (left column) and extended sources (right 
column). 



Dipole Fit Minimum LI Norm, max curr. 

Figure 10: MEG dipole fit result (left) and LI norm 
reconstruction with a maximum current density 
constraint (right). 




